Chiral dynamics and the growth of the nucleon's gluonic transverse size at small x 
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We study the distribution of gluons in transverse space in the nucleon at moderately small a; (~ 
10~^). At large transverse distances (impact parameters) the gluon density is generated by the "pion 
cloud" of the nucleon, and can be calculated in terms of the gluon density in the pion. We investigate 
the large-distance behavior in two different approaches to chiral dynamics: phenomenological soft- 
pion exchange, ti) the large- A*';, picture of the nucleon as a classical soliton of the pion field, which 
corresponds to degenerate N and A states. The large-distance contributions from the "pion cloud" 
cause a ~ 20% increase in the overall transverse size of the nucleon if x drops significantly below 
Mtt/Mn ■ This is in qualitative agreement with the observed increase of the slope of the t-dependence 
of the J/tp photoproduction cross section at HERA compared to fixed-target energies. We argue that 
the glue in the pion cloud could be probed directly in hard electroproduction processes accompanied 
by "pion knockout", 7* + A —> 7 (or p, J/^) + n + N' , where the transverse momentum of the 
emitted pion is large while that of the outgoing nucleon is restricted to values of order Af^ . 



I. INTRODUCTION 

Hard scattering processes induced by (real or virtual) 
photons are an important source of information about the 
structure of the nucleon. The parton model description 
of such processes distinguishes a "longitudinal" direction, 
defined by the 3-momentum of the incoming photon, and 
the "transverse" plane perpendicular to it. The structure 
functions of inclusive deep-inelastic scattering are pro- 
portional to the parton densities, describing the distribu- 
tion of partons with respect to longitudinal momentum 
of the parent nucleon; they do not carry any information 
about the distribution of partons in the transverse plane. 
In this sense, inclusive deep-inelastic scattering provides 
us with a 1-dimensional image of the nucleon. Much 
more detailed information can be obtained from exclu- 
sive processes, in which one measures amplitudes with 
a non-zero momentum transfer between the initial and 
final nucleon, A = p' — p ^ 0, with t = < 0. These 
include hard electroproduction processes such as deeply- 
virtual Compton scattering and meson production, or 
diffractive photoproduction of heavy quarkonia {J/ip, T), 
which probes the gluon distribution in the nucleon. Such 
measurements can give information also about the spa- 
tial distribution of partons in the transverse plane, thus 
providing us with a 3-dimensional image of the nucleon. 

A precise formulation of the notion of a spatial dis- 
tribution of partons in the transverse plane is possible 
within the formalism of generalized parton distributions 
(GPD's), which parametrize the non- forward matrix el- 
ements {p' ^ p) of twist-2 QCD light-ray operators be- 
tween nucleon states P, 0, l3„ ^^J . Of particular interest 
is the "diagonal" limit of zero longitudinal component 
of the momentum transfer, A+ = 0. The GPD's in 
this case depend (in addition to the partonic variable, 
x) only on the transverse component of the momentum 
transfer, A^, with — A5_ = t, and the nucleon helicity- 
conserving ones reduce to the usual polarized and un- 
polarized parton densities in the limit i — > (for this 
reason the A+ = GPD's are also referred to as "non- 



forward parton densities" ^]). The Fourier transform of 
these functions with respect to then defines func- 
tions of a 2-dimensional coordinate variable, b, which 
can be interpreted as the spatial distributions of par- 
tons with longitudinal momentum fraction x in the trans- 
verse plane ( "impact parameter-dependent parton distri- 
butions") Q- The integral of these b-dependent distri- 
butions reproduce the total densities of partons for given 
X, and they can be shown to satisfy positivity conditions 
locally in b This "mixed" momentum and coordinate 
representation provides a vivid 3-dimensional (more pre- 
cisely, 2-dimensional) picture of the partonic structure 
of the nucleon, which naturally lends itself to the visual- 
ization of polarization effects, including orbital angular 
momentum |^. On the phenomenological side, it has 
been speculated that this representation could provide a 
simple explanation of the transverse single spin asymme- 
tries observed in semi-inclusive deep-inelastic scattering 

s 

It should be noted that the "diagonal" GPD's (A+ = 
0) are not directly observable, since hard exclusive pro- 
cesses generally have skewed kinematics. However, in a 
number of important cases the kinematics is effectively 
not far from the diagonal case. In particular, this is the 
case for J/ip photoproduction, which probes the gener- 
alized gluon distribution in the nucleon; the ratio of the 
longitudinal momentum fractions of the two gluons is of 
the order Xi/x2 ^ 2 j9|. Furthermore, for T production 
[lOj , and generally for electroproduction at large and 
sufficiently small x j^, the gluon GPD's are dominated 
by evolution down from significantly larger values of xi 
and X2. Since the difference xi — X2 is not changed by 
the evolution, this implies that one is actually probing 
the ^-dependence of the nearly diagonal distributions at 
larger xi and X2 (the ^-dependence is not affected by the 
evolution) . 

An interesting aspect of the impact parameter repre- 
sentation of parton distributions, which has been less em- 
phasized so far, is that it provides a useful new frame- 
work for developing dynamical models for (generalized) 
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parton distributions. The possibility to separate contri- 
butions from different transverse distance scales allows 
one to introduce the notion of a "size" of the nucleon, 
well-known from elastic form factors (charge radii), in 
the discussion of generalized parton distributions. Fur- 
thermore, one can incorporate the fundamental fact that 
the long-distance behavior of strong interactions is gov- 
erned by chiral dynamics. Thus, one would expect that, 
in a certain range of a;, the parton distributions at trans- 
verse distances of order 1 /M^ can be attributed to the 
"pion cloud" of the nucleon, and can be computed from 
first principles in terms of the known parton distributions 
in the pion. This would provide model-independent con- 
straints for the distributions of partons in the transverse 
plane. Moreover, the partons in the pion cloud may lead 
directly to observable effects in hard processes sensitive 
to large transverse distances. 

In this paper we study the impact parameter depen- 
dence of the nucleon's gluon density from a phenomeno- 
logical perspective. We concentrate on the region of mod- 
erately small X (~ 10~^), for which the i-dependence 
of the gluon GPD (the so-called "two-gluon form fac- 
tor" of the nucleon) has been extracted from a number 
of electro-/photoproduction experiments, see for 
a compilation. The investigation consists of three parts. 
First, we derive the asymptotic behavior of the impact 
parameter-dependent distribution at large 6 = |b|. In 
this limit the gluon distribution is dominated by soft- 
pion exchange contributions and calculable in terms of 
the gluon density in the pion. The asymptotic behavior 
at large h can be stated in the form of a DGLAP-type 
convolution of the gluon density in the pion with a 5- 
dependent distribution of pions in the nucleon, which is 
completely determined by chiral dynamics and exhibits a 
Yukawa-like falloff at large h. Second, we investigate the 
contribution of large transverse distances, 6 ^ I/Mtt, to 
the average transverse size of the nucleon, (6^). The lat- 
ter is proportional to the slope of the i-dependent gluon 
GPD at t = 0, and thus directly observable. We find that 
the large-distance contributions lead to a finite increase 
in the transverse size of the nucleon if x drops signifi- 
cantly below ~ Mj^/Mm- Using a simple two-component 
picture, we estimate the magnitude of the increase at 
roughly 20% of the total transverse size of the nucleon. 
Such an increase is indeed observed when comparing the 
slope of the i-dependence of the J /ip photoproduction 
cross section at HERA with that at fixed-target ener- 
gies Third, we ask how the glue in the pion cloud 
could be probed directly in hard exclusive processes. A 
promising candidate is hard electroproduction of photons 
(or p, J /ip mesons) accompanied by the "knockout" of an 
additional pion, 7* -I- TV 7 (or p, J/^) -I- tt -t- TV', 
such that the transverse momentum of the pion is large 
while that of the outgoing nucleon is restricted to values 
of order A/^r. 

It is worth noting that there is an interesting connec- 
tion between the increase of the transverse size of the 
nucleon due to the pion cloud described here, and the 



so-called shrinkage of the diffractive slope in soft physics. 
In soft physics the shrinkage is usually interpreted as due 
to Gribov diffusion of the impact parameters in the par- 
tonic ladder, see Ref. 13] for a pedagogical discussion. 
The phenomenon we discuss here is essentially related 
to the diffusion in the first rung of the ladder, which is 
delayed as compared to soft physics by selection of rela- 
tively large-a; gluons which are missing in the pion cloud. 

A crucial element in our approach is the impact- 
parameter dependent "parton distribution" of pions in 
the nucleon at large transverse distances, h > I/A/tt- We 
show that this concept is meaningful only for pion mo- 
mentum fractions y < Mt^IMj^^ where the pion virtual- 
ity is of order Af^. We calculate this distribution in two 
different approaches to chiral dynamics — soft-pion ex- 
change with a phenomenological pion-nucleon coupling, 
and the large-A'^c picture of the nucleon as a classical 
soliton of the pion field. We explicitly demonstrate the 
equivalence of the two approaches. The key to this is the 
realization that in the large-iVc limit the nucleon and the 
Delta resonance are degenerate, so the large- A'c calcu- 
lation based on the classical soliton corresponds to soft- 
pion exchange with iV and A intermediate states included 
on the same footing. We also comment in general on the 
relation of pion cloud contributions at large impact pa- 
rameters to the 1 /A^c~expansion of parton distributions 
[1^ . Our use here of the large-TVc limit is completely gen- 
eral, and not related to any particular dynamical model 
used to generate stable soliton solutions. Nevertheless, 
the results obtained here provide useful constraints for 
calculation of parton distributions within specific l arge - 
models, such as the chiral quark-soliton model 

Superficially, our approach bears some resemblance to 
the well-known pion cloud model used to explain the ob- 
served flavor asymmetry of the sea quark distributions in 
the proton, d(x) > u{x) 0, see e.g. Refs. 17] for recent 
reviews. There is, however, an important conceptual dif- 
ference in that we invoke the notion of pion exchange only 
at large transverse distances, b ^ l/Af^, where the con- 
tributions from nN configurations in the nucleon wave 
function are distinct from those of average configurations 
and thus physically meaningful. This way of defining the 
"pion cloud" contributions naturally avoids a number of 
problems which have plagued the traditional pion cloud 
model, see e.g. the discussion in Ref.|l8l]. In our ap- 
proach the pions are almost on-shell (the virtuality is of 
order Af^), we do not need cutoffs at the pion-nucleon 
vertices, and the results are independent of the form of 
pion-nucleon coupling (pseudoscalar or axial). 

In the present investigation we concentrate on the un- 
polarized gluon distribution, the i-dependence of which 
is well-known in the regions of x covered by J/4' pho- 
toproduction at fixed-target and collider energies p^ . 
The approach developed here can easily be extended to 
the quark distributions (singlet and non-singlet), as well 
as to polarized quark and gluon distributions. It can 
also be applied to the diagonal (zero skewedness) limit 
of GPD's with a helicity fiip between the nucleon states, 
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which have no correspondence in the parton distributions 
of inclusive deep-inelastic scattering. 

The study of the transverse size of the gluon distri- 
bution in the nucleon is important also for a number 
of other problems, not directly related to exclusive pro- 
cesses. Knowledge of the transverse size of the gluon 
distribution at moderate x is required for modeling the 
initial conditions for non-linear QCD evolution equations 
descri bing the saturation of parton densities at very small 
X 0, 120I I2H . The transverse size of the hard compo- 
nent of the nucleon wave function is also an important 
parameter in modeling hadron-hadron collisions at LHC 
energies [2^ . 

This paper is organized as follows. In Section^ we re- 
view the basic properties of the impact parameter repre- 
sentation of the gluon density, starting from its covariant 
definition as the Fourier transform of the i-dependent 
non- forward gluon density. In Section IlIII we derive the 
asymptotic behavior of the gluon density in the nucleon 
at large impact parameters as due to soft-pion exchange, 
and discuss its dependence on x. In Section llVl we address 
the same problem from the point of view of the large- A^^ 
limit, where the nucleon is described as a classical soliton 
of the pion field, and demonstrate the equivalence of this 
picture to soft-pion exchange including Delta resonance 
contributions. In Section we show that the pion cloud 
contribution at 6 ^ produces a finite increase of the 

overall transverse size of the nucleon, {b'^), if x drops sig- 
nificantly below Mt^/Mn, and estimate the magnitude 
of the increase in a simple two-component picture. In 
Section IVII we explore the possibility of measuring the 
glue in the pion cloud directly by way of electroproduc- 
tion of photons or vector mesons at small t accompanied 
by large-angle pion production ("pion knockout"). Our 
conclusions and possible generalizations of the approach 
proposed here are summarized in Section IVIII 



II. 



IMPACT PARAMETER REPRESENTATION 
OF THE GLUON DENSITY 



The matrix element of the twist-2 QCD gluon opera- 
tor between nucleon states of different momenta exhibits 
two Dirac structures (nucleon helicity-conserving and - 
flipping), and is parametrized as 



dX 
2^ 



X G^^(-An/2) G^^(An/2) \N{P - A/2)) 
= xHgix,t)u'hu + ^^'/j^:^^ uV,..A^n-u. (1) 



2Mn 

Similar definitions apply to the matrix elements of the 
twist-2 quark operators, see Ref.0- Here, n is a light- 
like four-vector, = 0, whose normalization is arbi- 
trary, and G^^ denotes the gluon field strength (here and 
in the following, a gauge link between the fields is implied 



but will not be written). The incoming and outgoing nu- 
cleon four-momenta are expressed in terms of the average 
momentum, P, and the momentum transfer, A. The nu- 
cleon spinors, u = u{P — A/2) and u' = u{P + A/2) (we 
suppress the helicity labels for brevity) are normalized 
according to uu = u'u' = 2Mj^, and our conventions are 
n = n^^^fi etc., and = (l/2)[7^,7^]. The functions 
Hg and Eg in Eq.JQ) depend on the partonic variable, x, 
and the invariant momentum transfer, t = IS? < 0, and 
are referred to as the non-forward gluon densities of the 
nucleon. It is implied here that 



(An) 
JPn) 



0; 



(2) 



in the general case the functions Hg and Eg would de- 
pend also on the "skewedness" parameter, ^. Finally, the 
functions Hg and Eg depend implicitly also on the nor- 
malization point of the QCD operator; this dependence is 
governed by evolution equations. In what follows we have 
in mind a typical scale for hard processes, say ~ 4 GeV^, 
which will not be exhibited explicitly. 

The non-forward gluon densities are defined on the in- 
terval — 1 < x < 1. Due to G-invariance they are odd 
functions of x, 



Hg{~x,t) = -Hg{x,t), 



(3) 



and similarly for Eg{x,t). At t = the function Hg 
for a; > reduces to the usual gluon distribution in the 
nucleon (as defined e.g. in Ref . j23l| ) 



Hg{x,t = 0) 



for X > 0. 



(4) 



The n'th moments in x of the functions Hg(x,t) and 
Eg{x,t) define the form factors of the local twist-2 spin- 
n gluon operators, which arise from the Taylor expansion 
of the non-local light-ray operator in Eq. |^ in powers of 
the separation, A. In particular, the second moments 
are related to the form factors of the gluonic part of the 
QCD energy-momentum tensor |2|. Note that with our 
conventions 



- / dxxHg{x,t — 0) — f 
2 J-i Jo 



dxxg{x) — Ag, (5) 



where Ag parametrizes the traceless part of the expec- 
tation value of the gluonic part of the QCD energy- 
momentum tensor in the nucleon, 



{N{P)\ Gf^(0)G^^(0) \N{P)) 



2Ag P^Pp 
trace terms, (6) 



and can be interpreted as the total momentum fraction 
carried by gluons in the parton model. 

For the discussion of the t-dependence of the non- 
forward gluon density it is convenient to define the so- 
called two-gluon form factor of the nucleon 0, 0| 



r(x,t) 



Hg{x,t=0) 



[x > 0), 



(7) 
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which is normahzed to T{x,t = 0) = 1. This function 
can be interpreted as the form factor describing the dis- 
tribution of partons with given longitudinal momentum 
fraction, x, and can be directly compared with the well- 
known elastic form factors of the nucleon (electromag- 
netic, axial). 

The non-forward gluon densities are invariant func- 
tions, not associated with any particular reference frame. 
However, as with the invariant elastic form factors, it is 
possible to give an interpretation of these functions as 
Fourier transforms of certain "charge densities" in a spe- 
cial frame in which the momentum transfer has only spa- 
tial components (6]. Choosing n to define the light-cone 
"plus" direction, ^ (1, 0, 0, -1), 



P- 



(nP) = P" 

pO _ p3 



pi 



(8) 
(9) 



one goes to a frame in which Pj^ = 0. The condition 
(riA) = implies A+ ~ 0, and the kinematical constraint 
(PA) — requires that also A~ = 0, i.e., the momentum 
transfer has only a transverse component, A_l, with 

Ai = ~t. (10) 

In this frame the non-forward density can be thought of 
as a function of x and the transverse vector Aj.. The de- 
pendence on Aj_ can be represented as a Fourier integral 
over a transverse coordinate variable ("impact parame- 
ter"), b: 



if,(x,-A2 



(11) 



where b = |b|. The limiting relation Q for the non- 
forward gluon density at t = then implies that 



(fbg{x,b) = g{x). 



(12) 



Thus, one may think of the function g{x^ b) as describing 
the distribution of partons of given longitudinal momen- 
tum fraction x in the transverse plane, in the sense that 
the integral over the transverse plane gives the usual par- 
ton densities, depending only on x. In Ref.|^ the name 
"impact parameter-dependent parton distribution" was 
proposed for these functions. 

An important feature of the impact parameter-depen- 
dent parton distributions is that they satisfy positivity 
conditions locally in b 7J. In particular, it makes sense 
to compute averages over the transverse plane, with the 
function g{x,b) acting as a positive definite weight. For 
instance, one can define the (cc-dependent) average trans- 
verse size squared of the nucleon as 



. ^ jd%g{x,b) b^ 
~ Jd%g{x,b) 



(13) 



the integral in the denominator is just the total gluon 
density, 5(2;). This average can also be expressed as the 



i-derivative of the non-forward gluon density, viz. the 
two-gluon form factor, Eq.©, at t — 0. Differentiating 
Ea. pifl twice with respect to the vector Aj^ at A^ = 0, 
and making use of the relation (|10|l . one obtains 



Hg{x,0) 



dHg{x,t) 
dt 



= 4 



dT{x,t) 
dt 



(14) 



The factor 4 here represents two times the dimension 
of the transverse plane; it replaces the factor of 6 in the 
well-known relation between the slope of the elastic form 
factors and the "three-dimensional" charge radii. 



III. LARGE TRANSVERSE DISTANCES AND 
THE PION CLOUD 

We now study the behavior of the nucleon's gluon dis- 
tribution at large impact parameters from a phenomeno- 
logical point of view. For understanding the physical 
mechanism governing the large-6 asymptotics one may 
start with the moments in x of the i-dependent nonfor- 
ward gluon densities, i.e., the nucleon form factors of 
the local twist-2 spin-n gluon operators, which are func- 
tions only of the invariant momentum transfer, t. Their 
Fourier transforms with respect to Aj_ (with t — — A5_) 
define the moments of the impact parameter-dependent 
gluon distributions. The asymptotic behavior of the lat- 
ter at large b is governed by the singularities in t of the 
i-dependent moments. Since the local twist-2 gluon op- 
erators are G-parity-even, the singularity closest tot = 
is the cut at i > 4M| coming from two-pion exchange in 
the t-channel. The threshold behavior of the imaginary 
part sX t ^ 4M^ -I- dominates the large-6 behavior of 
the moments of the 6-dependent gluon distribution. 

Having noted this, we may study the pion exchange 
contributions to the non-forward gluon densities directly 
for the x-dependent functions, i.e., for the matrix ele- 
ment of the non-local gluon operator of Eq. iQJ , see Fig. 
[37I I. The dashed lines indicate the pion, the solid line the 
nucleon (later we shall consider also graphs with A in- 
termediate states), and a phenomenological pion-nucleon 
coupling is assumed at the vertices. Before specifying 
how to interpret and evaluate the graphs of Fig. ^ let 
us emphasize that we shall use them only as a means to 
derive the large-6 asymptotics of the impact parameter- 
dependent gluon distributions in the nucleon (or, equiv- 
alently, the t 4M^ -I- threshold behavior of the imag- 
inary part of the t-dependent non-forward gluon densi- 
ties), which is governed by exchange of soft pions, much 
like the nucleon-nucleon interaction at large distances 
(Yukawa potential). We shall not attempt to make sense 
of the resulting Feynman integrals in regions where the 
pion virtuality is >• M^. 

The blob in the graphs of Fig. ^ represents the non- 
forward gluon density in the pion. In the physical region 
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-An/2 An/2 



-Xtl/I XtlI'I the electromagnetic form factor, 




TT 



TT 



+ 



N 



N{A) 



N 




FIG. 1: The pion exchange contribution to the non-forward 
gluon densities of the nucleon, Eq.Q, determining the lead- 
ing asymptotic behavior of the impact parameter-dependent 
gluon densities at large b. 



t < this function parametrizes the on-shell matrix ele- 
ment of the QCD gluon operator between physical pion 
states, 

n'^n^ {n'^ik + A/2)| G^T(-An/2) G^^(An/2) 
X |7r''(fc - A/2)) 

^ S"'' {knf f dz e-'''^^'''''' zHg,^{z,t), (15) 



where again (An) — 0. The non-forward gluon den- 
sity in the pion Hg^T^{z,t) is an analytic function of t 
which can be continued to the unphysical region. In 
the soft~pion exchange contribution to the distribution in 
the nucleon this function enters at the two-pion thresh- 
old, t — 4M^ > 0, in the unphysical region. By cross- 
ing symmetry, it could be related to the two-pion gluon 
distribution amplitude (i.e., the matrix element of the 
same gluon operator between the vacuum and a two-pion 
state) in the physical region for two-pion production; for 
a discussion of the crossing relations see Befs.p^j. In 
the following we shall neglect the ^-dependence of the 
non-forward gluon density in the pion, i.e., we replace 
the nonforward by the usual {t — Q) gluon density in the 
pion, g^{z): 



Hg^^{z,t) 



Hg^^{z,t^O) 



g^{z), z > 0, 
-g^{~z), z <0. 

. (16) 

To justify this approximation, we note that, while the 
exact i-dependence of the non-forward gluon density in 
the pion is unknown, its characteristic scale can be safely 
assumed to be much harder than the 4M^ involved in 
the extrapolation from t — to the two-pion threshold. 
A reasonable guess is to suppose that the distribution 
of gluons in the pion follows that of the quarks, which 
would imply that the two-gluon form factor of the pion 
[defined in analogy to Eq.jTI)] can be approximated by 



1 - t/m'^ 



2ff 



'2s 



Mi 



0.5 GeV^ 



(17) 

(we assume the form factor to be independent of the 
gluon momentum fraction z, in the region of interest). 
The analogy with the nucleon, where the two-gluon form 
factor has been extracted from the i-dependence of J/ip 
photoproduction data, suggests that its scale may be 
even harder, m^g ~ 1 GeV^ 12]. Either way, one con- 
cludes that the continuation from < = to 4A/^ in the 
unphysical region should have a small effect (~ 10%). 

Another way of motivating the approximation Eq. Hlti|) 
is to say that we shall assume that the intrinsic trans- 
verse size of the pion is small compared to the transverse 
distances, 6, at which we want to study the gluon distri- 
bution in the nucleon. This is certainly justified as long 
as we are interested in the asymptotic behavior in the 
strict b —^ oo limit. Later, when we use the asymptotic 
expression at finite b to estimate the pion cloud contri- 
bution to the nucleon's transverse size, we shall expect 
corrections from the intrinsic transverse size of the pion, 
see the discussion in Section Ivl 

The above approximation can be stated in an alter- 
native way, by defining an operator in the interpolating 
pion field, which represents the QCD gluon light-ray op- 
erator in the effective chiral theory. One easily sees that 
Eas. H15|l and (|16|l amount to matching the QCD gluon 
operator with a pionic light-ray operator as 

n^n^* G^^^Xn/2)G^p{Xn/2) 

-n°'n'^ I dzzg^{z) 



Ti°'(-z\n/2) da d0TT°'{z\n/2) + {z ^ -z) 



-\- total derivatives. 



(18) 



where 5= (l/2)(9 — d)- We have omitted here terms 
involving total derivatives contracted with the light-like 
vector n, which contribute only for non-zero longitudinal 
momentum transfer, c/. Eq.Q. Computing the matrix 
element of the pionic operator (|18|l between pion states 
and comparing with the general parametrization Ea. p5|l . 
one obtains i?g^^(z,i) = g-n{z) for z > 0, in agreement 
with Ea. (|16(l . The approximation of neglecting the t- 
dependence of the gluon distribution in the pion is now 
encoded in the fact that the pionic operator does not 
contain any contracted total derivatives (or, equivalently, 
no transverse total derivatives) , which would give rise to 
powers of A^^ = —t in the matrix element j^. The 
pionic operator is defined also off-shell, to the extent that 
the above approximations are valid. In particular, with 
the operator representation of the gluon density in the 
pion, we are allowed to interpret the graphs of Fig. ^ as 
Feynman graphs, which can be evaluated using standard 
methods. 
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A comment is in order concerning chiral invariance. 
Strictly speaking, the QCD gluon operator should be 
matched with a chirally invariant pionic operator, con- 
structed from fields that transform according to a non- 
linear realization of the chiral group, SU{2)l x SU{2)b.- 
The bilinear pion operator of Eq. (|18|l is not chirally in- 
variant; however, it corresponds to the small-field limit 
of such a chirally invariant operator. This simplification 
is justified in our case, since we shall only evaluate the 
pionic operator in the region of large impact parameters 
h ^ l/A^TT, where the pion field is exponentially small 
in h. Our results for the asymptotic behavior for large 
h below would be the same if we started from a chirally 
invariant, non-linear operator. 

Having at hands the representation of the QCD twist- 
2 gluon operator as a pionic operator, Ea. (|18|) . it is 
straightforward to compute the pion exchange contribu- 
tion to the non-forward gluon densities in the nucleon. 
Substituting the pionic representation of the gluon oper- 
ator in the matrix element of Eq.J^I, and changing the 
integration variable to y = x/ z, one obtains (here a; > 0) 

Hg{x, t)cioud = / — grr{x/y) H^{y,t), (19) 

Jx V 



Eg (x, t)cloud 



1 dy 



g^{x/y)E^{y,t). (20) 



The functions and i?^ parametrize the nucleon ma- 
trix element of the pion light-ray operator, in complete 
analogy to the gluon distributions of Eq. ^2^1 ; 



(7V(P + A/2)| 



X 7r''(-An/2)7r'^(An/2) |A^(P- A/2)) 
= H^{y,t)u'hu + E^!^u'a^,^^^n'^u. (21) 



2Mjv 

They can be interpreted as the isoscalar non-forward 
"parton densities" of pions in the nucleon, with y the lon- 
gitudinal momentum fraction of the pion. We emphasize 
that this notion is meaningful only in the region in which 
the pion virtualities are of order ; we shall see below 
that this is indeed the case for y < M-^/Mn. In the im- 
pact parameter representation [cf. Ea. (|ll|l ]. Eas. p9|l and 
(Pn)l take the form (here x > 0) 



loud 



gA^/y) Uiy.b), (22) 

X V 



and similarly for the helicity flip-distributions. The b- 
dependent distributions of pions in the nucleon are de- 
fined in analogy to Ea. (|ll|l . 



F,(y,^Ai) = j d^ 



f7r{y,b), 



(23) 



and similarly for the helicity-flip distribution ET^{y,t). 
Again, these distributions are meaningful only at dis- 
tances b ^ I/Mtt, where the pion virtualities are of order 



M^, as will be shown below. We refer to Ea. H22() as the 
"pion cloud contribution" to the gluon density in the nu- 
cleon. Note that the relation between the distribution of 
pions and of giuons is local in impact parameter space; 
this is a result of our approximation of neglecting the 
intrinsic transverse size of the pion. 

The distribution of pions in the nucleon at impact pa- 
rameters b > I/A/tt is determined by chiral dynamics 
and can be computed from first principles, regarding the 
pion and nucleon as pointlike particles and using the phe- 
nomenological ttNN coupling. For the moment we con- 
sider only the graphs in Fig.nwith intermediate nucleon 
states, which dominate in the strict b —^ oo limit; contri- 
butions with intermediate A states will be discussed in 
Section Hvl in connection with the large-iVc limit. The 
essential steps of the calculation are presented in Ap- 
pendix^ One computes the t-dependent distribution of 
pions through the matrix element of pion light-ray oper- 
ator in Ea. (|21|l . and then studies the asymptotic behavior 
of its Fourier transform at large b. This can be done ei- 
ther by writing the Feynman integral in an appropriate 
parameter representation f Appendix lA l|l . or by directly 
computing the imaginary part of the i-dependent distri- 
bution using the Cutkosky rules (Appendix lA 2|l . The 
distributions at large b are finite without cutoffs at the 
pion-nucleon vertices, and independent of the type of 
pion-nucleon vertex (pseudoscalar or axial vector) em- 
ployed in the calculation. The result for the asymptotic 
behavior at 6 — s- cx) is 



f-r^iy^b) 



igl 



167r 



NN y_ 
y 



2y^Ml 



Ml 



-KNb 



KNb 



UN{y,b). 



(24) 



To indicate the fact that this part of the distribution of 
pions is obtained from intermediate nucleon states only, 
we denote it by /ttAt, in order to distinguish it from the 
contribution from intermediate A states to be introduced 
below; a similar notation we adopt for the t-dependent 
distributions of Ea. (|21|) . Here and in the following. 



y = 



(25) 



At large 6 the distribution of pions drops exponentially, 
with the "decay constant" 



KM 



y 



Ml 



1/2 



(26) 



Note the similarity with the Yukawa potential in the NN 
interaction. However, in our case the decay constant de- 
pends on the longitudinal momentum fraction of the pi- 
ons in the nucleon, y. This j/-dependence, which arises in 
an elementary way from the singularities of the Feynman 
graphs in Fig. (see Appendix . plays a crucial role 
for the effects discussed in the following. 

The inverse decay constant, 1/kjv, is a measure of the 
transverse size of the distribution of pions in the nucleon. 
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FIG. 2: The inverse "decay constant", 1/kjv, of the exponen- 
tial tail of the distribution of pions in the nucleon at large 
impact parameters, Egs. 12411 and 1261 . as a function of the 
pion momentum fraction, y. This quantity determines the 
transverse size of the pion cloud for given y. The dashed line 
marks the value y = M^/Mjv — 0.14. The limiting value of 
1/kjv for y ^ is l/(2Af^) = 0.73 fm. Note that the notion of 
pion cloud contributions to the nucleon's parton distributions 
is physically meaningful only for pion momentum fractions 
y < Mtt/Mn, where the size of the pion cloud is significantly 
larger than those of average configurations in the nucleon. 



Its dependence on y exhibits two different regimes, see 
Ea. (PS|l . For y < M^^/Mm the size becomes independent 
of the nucleon mass and approaches the value 1 / (2M^) at 
y — Q. In this region the transverse size of the pion distri- 
bution is determined by the pion mass only (it would be 
singular in the chiral limit, Af^ 0). For y > M.^/Mn 
the transverse size is of the order of the inverse nucleon 
mass, 1/Mn- It is a slowly decreasing function of y, van- 
ishing linearly in 1 — y for y — > 1, in agreement with gen- 
eral expectations. Fig. [3 shows the transverse size 1/kjv 
as a function of y obtained with the physical values of the 
pion mass; the dashed line marks the value y = Mt^/Mj^, 
the boundary between the two regimes. The above im- 
plies that the very notion of pion cloud contributions to 
the nucleon's parton distributions is meaningful only for 
pion momentum fractions y < M^/M^r, where the size of 
the pion cloud is significantly larger than that of average 
configurations in the nucleon. 

The change of the spatial size of the pion distribu- 
tion with y is matched by a corresponding change of the 
pion virtualities. From Eas. HAQI) . (|A10|1 and (|Alip in 
Appendix^ one sees that the minimum value of the vir- 
tuality (taken with positive sign) of the pions in the dis- 
tribution Ea. H21|l . i.e., the virtuality of the pion lines in 
the graphs of Fig. ^ is 



N 



(27) 



see also the discussion in Ref . [isj . where this quantity is 
denoted by — imin- For y of order unity the pion virtuality 





1/M, 



N 



N 



1/M. 



FIG. 3: Visualization of the characteristic transverse size of 
the pion cloud contribution to the gluon density in the nu- 
cleon, as defined by Eqs.l|221l, IMl and x > Mt,/Mn. 
The transverse size of the pion cloud contribution (light 
shaded disc) is of order 1/Mjv, comparable to that of aver- 
age configurations in the nucleon wave function producing the 
"bulk" of the gluon distribution (dark shaded disc). In this 
regime the notion of pion cloud contributions is not meaning- 
ful. X <^ Mtt/Mm'. The size of the pion cloud contribution 
grows to I/Mtt. Now the pion cloud contribution is distinct 
from the "bulk" of the gluon distribution. The result in an 
increase of the overall transverse size of the nucleon (see Sec- 
tion |V|l. 



is of the order M^, and the notion of pion exchange 
becomes meaningless. It is only for values y < Mt^/Mm 
that the pion virtuality reduces to values of the order 
A/^, where the pion can be regarded as soft. As we see 
from the above discussion, this is the region where the 
transverse size of the pion distribution, 1/kjv, is of the 
order l/M^r. 

The transverse size of the distribution of pions in the 
nucleon determines that of the pion cloud contribution 
to the gluon density, Ea. H22|) . The convolution integral 
in Ea. H22() runs over values y > x. For x > Mt^/Mm 
one integrates only over the region where the size of the 
distribution of pions in the nucleon is ~ l/M^r, and the 
size of the pion cloud contribution to the gluon density 
will be accordingly. For x < Mj^/Mn, however, the in- 
tegral extends also into the region y < M^^/Mm where 
the size of the distribution of pions is ^ l/Af^r, so the 
size of the pion cloud contribution to the gluon density 
grows to ~ Note that for the gluon distribution at 

distances 6^1 /A/tt to become sizable one really needs x 
significantly smaller than M^^/Miq , so that in addition to 
y also the momentum fraction of the gluon in the pion, 
z — x/y, can become small. To summarize, the qualita- 
tive change in the transverse size of the distribution of 
pions between y > M^^/Mn and y < Mj,/Mn [cf. Fig.|21) 
directly translates into a corresponding change of the size 
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of the pion cloud contribution to the gluon density in 
the nucleon. This behavior is schematicaUy illustrated in 
Fig. 13 where the transverse size of the gluon density is 
represented by a disc in the transverse plane. The dark 
shaded disc indicates contributions from average configu- 
rations in the nucleon wave function to the gluon density, 
whose transverse size is of order l/A/jy. The light shaded 
disc indicates the pion cloud contributions, as defined by 
Eq.(I12l. For X > Af^/A^jv the pion cloud contributions 
are "hidden" in the average configurations. In this re- 
gion of X the asymptotic expression (|22|1 is only of sym- 
bolic significance, since the chiral contributions to the 
gluon density are not special in any sense. It is only for 
X <^ Mt^/Mm that the pions are free to propagate over 
transverse distances of order and the pion cloud 

contributions become a distinct feature of the gluon den- 
sity in the nucleon. As we have seen above, this is also 
the region where the pion virtuality is of order M"^, and 
our approximations are justified. The existence of these 
two different regimes has important consequences for the 
overall average transverse size squared of the nucleon, as 
will be shown in Section Ivl 

To conclude this discussion, we return to the question 
of the role of the intrinsic transverse size of the pion. 
The i-dependence of the non-forward gluon density in 
the pion given by Ea. (|17|l (assumed to be independent 
of the gluon momentum fraction, z, in the region of in- 
terest) would imply an exponential decay of the impact 
parameter-dependent gluon distribution in the pion at 
large impact parameters (with respect to the center of 
pion) of the form 



-m2gb' 



/32tt y/m2gb' 



(28) 



The characteristic size of this distribution, 



'25 



0.2...0.3fm, is smaller than the size of the distribu- 
tion of pions in the nucleon in the region y <C Mt^/Mn, 
l/(2Af^) — 0.73 fm, by a factor of ^ 3. This ratio gives 
an estimate of the accuracy of our "two-scale picture" , 
cf. Fig. 13 Note that when computing the average trans- 
verse size squared of the nucleon in Section^the relevant 
ratio is that of the transverse areas covered by the distri- 
butions, so we can expect a reasonable accuracy in this 
case. 



IV. DISTRIBUTION OF PIONS IN THE 
LARGE-TVc LIMIT 

So far we have studied the asymptotic behavior of the 
gluon distribution at large impact parameters invoking 
the phenomenological notion of soft pion exchange. It is 
interesting to investigate the large-6 behavior also in a 
different but related approach to chiral dynamics, namely 
the large- A^c limit of QCD, where the nucleon is de- 
scribed as a classical soliton of the effective chiral the- 
ory. It turns out that in this approach one can relate the 



distribution of pions in the nucleon at large b directly to 
the large-distance "tail" of the classical pion field in the 
soliton rest frame, which allows for a simple geometric 
interpretation. The functional form and the magnitude 
of this "tail" are completely determined by chiral sym- 
metry, and thus independent of the particular dynamical 
model used to generate stable soliton solutions. This al- 
lows us to conduct the following discussion at a general 
level, assuming only a "generic" chiral soliton picture of 
the nucleon. 

We begin by studying the i-dependent distribution of 
pions in the nucleon, Ea. (|21|l . in the large-TVc limit. Fol- 
lowing the standard procedure for form factors, the 1/Nc 
expansion is performed in the Breit frame, in which the 
three-momenta of the incoming and outgoing nucleons 
are 0(1) in N^, and thus suppressed compared to the 
energies, which are 0{Nc) [the nucleon mass, M^, is 
0{Nc)]. For the average momentum P and the momen- 
tum transfer A [cf. Eq.l^] this implies 



pO 

AO 
A* 



Mn 

0, 

0, 

0(1) 



o(i), 



(29) 
(30) 
(31) 
(32) 



Choosing the spatial component of the light -like vector 
n along the 3-direction, — (1,0,0,-1), the condition 
that (nA) = requires A^ = 0, i.e., A^^ = = 0(1). 
Since the impact parameter b is Fourier-conjugate to 
Aj^, this implies that we are studying the transverse spa- 
tial distributions at distances b = 0(1) in 1/Nc. This is 
natural in view of the fact that typical nucleon radii (elec- 
tromagnetic, axial) are 0(1) in the large-iVc limit. Con- 
cerning the pion longitudinal momentum fraction, since 
M-^/Mn — 1/iVc we suppose that we compute the distri- 
bution of pions in the nucleon for "average" values of the 
pion momentum fraction, y ^ l/Nc- 

It is then straightforward to compute the matrix ele- 
ment in Ea. (|21|l in leading order of the 1/Nc expansion 
in a "generic" chiral soliton picture of the nucleon. One 
evaluates the pionic operator in the classical pion field 
of the soliton and projects on nucleon states of definite 
spin/isospin and momentum by integrating over collec- 
tive (iso~) rotations and translations of the classical soli- 
ton with appropriate collective wave functions 27]. Tak- 
ing the Fourier transform with respect to the momen- 
tum transfer A^ we obtain a simple expression for the 
isoscalar distribution of pions at large impact parame- 
ters, 

-OO ■^'^ J OO 

^ ""asymp (''-)''''asvmp(''+)i (^3) 



where 



rzr = b + =F 77 ^3 



(34) 



9 



trans V. 




longitud. 



FIG. 4: Schematic illustration of the calculation of the dis- 
tribution of pious in the nucleon at large impact parameters 
in the large-A^c limit, /7r(j/, 6)iargc-Arc' Eg. (1331 . Shown is a 
2-dimensional projection of 3-dimensional coordinate space 
in the soliton rest frame, containing the longitudinal and one 
transverse direction (say, the 3-1 plane). The classical pion 
field of the soliton is indicated by the circles; the direction 
of the field in isospin space at a given point is determined 
by the radius vector. The positions r_ and r+ at which the 
pion field is evaluated are indicated by crosses. Their separa- 
tion in the longitudinal direction is A. When integrating over 
the longitudinal coordinates A and uj, the points move along 
the line indicated in the figure, which passes the soliton at a 
transverse distance b from the center. 



are three-dimensional coordinate vectors in the frame 
where the soliton is at rest and centered at the origin, 
and 



asymp 



(r) 



'ig-nNN ( 1 



r 



(35) 



is the "Yukawa tail" of the classical soliton at large dis- 
tances, which is obtained from solving the linearized field 
equations for the pion field, see e.g. Ref . . A visual 
representation of the integral 1)33(1 defining the distribu- 
tion of pions is given in Fig. ^ At large b, which is where 
Ea. ((33|l applies, the integrals over the longitudinal coor- 
dinates A and u) can be computed in the saddle point 
approximation. For the leading asymptotic behavior we 
find 



/7r(2/, 6)largc-Afe — 



95L_ 
167r2 



where 



2 2/2^/2, + 



.1/2 



(36) 



(37) 



When comparing the large-iV^ result, Eas. (|^ and 
(|37|l . with the result of the phenomenological soft-pion 



exchange graphs, Eas. (|24|l and (|26() . we need to take into 
account that the large- A^c result implies y ^ <C 1. 

Indeed, we see that for y <ti 1 one has Kqo = kn, and the 
functional forms of Eqs.|(23J) and (|36|) coincide. However, 
we observe that the coefficient of the large- A^c result is 
larger than that of the pion loop calculation by a factor 
of 3: 

/7r(2/, &)largc-W, = 3/^Ar(y, 6) (y 1/A^c)- (38) 

This paradox is resolved by noting that in the large-A^c 
limit the nucleon and the Delta resonance are degenerate. 
A/a — Afjv ^ l/A'c. Thus, the pion loop calculation in this 
case should include A intermediate states on the same 
footing as nucleon ones, i.e., the relevant distribution of 
pions is now 



Uiy,b) = UNiy,b) + UA{y,b), 



(39) 



where f^rAiy, b) denotes the contribution from the graphs 
in Fig. n with A intermediate states. 

To verify the equivalence of the two approaches explic- 
itly, we compute the A contribution to the distribution 
of pions at large b in the approach of Section llTTl describ- 
ing the A resonance by an elementary Rarita-Schwinger 
field, and assuming a phenomenological ttA^A vertex. De- 
tails are presented in Appendix lA 31 In the general case. 
Ma > Mm and < y < 1 arbitrary, we obtain for the 
leading large-5 behavior 



UA{y,b) 



rj2 



Ml 



yM. 



N 



-(Mat- Ma)' + 



Ml 
2 



2y, 



y 

4Ml^ 
-KAb 



M% 



y 

Ml? 



KAb 



{Mi - M^) 



(40) 



where 



KA 



yHdj, ^ y{Ml-Mj,) 



Mi 



y 



y 



1/2 



(41) 

Note that for Ma > Mm the exponential decay of the 
A contribution is faster than that of the A^ contribution, 
c/. Eas. l|24|) and H26|l . In the large- A'^j, limit, however, 
where Mat, Ma ~ N^, Mn - Ma l/Nc, and M^ ~ 1, 
and for y ~ 1/A^c, we can neglect the y{M\ — M'j^)/y 
term in Eq.l^lJ, and thus have ka = f^N + 0{1/Nc), 
i.e., the two contributions exhibit the same exponential 
decay. Using the same relations we can simplify also the 
pre-exponential factors. Finally, taking into account that 
in the large- A^c limit 



zgTrNN, 



(42) 



10 



we see that the A contribution in this Hmit is exactly 
twice as large as the N contribution, Ea. (|24|l : 

UA{y,b) = 2UN{v,b) (large-iVj. (43) 

The sum of the N and A contributions, Ea. H39l) . in the 
large-A'c limit is thus exactly 3 times the N contribution 
alone, and coincides with the result of the calculation 
using the classical soliton field, Ea. (|36|l . 

To summarize, in the large- A^c limit the calculation of 
the distribution of pions from the "tail" of the classical 
soliton and from phenomenological soft-pion exchange 
graphs, including both N and A intermediate states, give 
the same result. Due to the degeneracy of the N and A 
the large-iVc limit leads to a different asymptotic behav- 
ior at large h compared to an "exact" phenomenological 
approach based on physical states: In the "exact" ap- 
proach the A contributions would be exponentially sup- 
pressed relative to the N ones and thus formally sub- 
leading, while in the large-A^c limit both come with the 
same exponential factor and contribute on the same foot- 
ing, with the A contribution even dominating due to the 
larger spin-isospin degeneracy. The fact that the A con- 
tribution is sizable in the large-TVc limit suggests that it 
is quantitatively important also in the "exact" approach 
and should be included, even though it is formally sub- 
leading. We shall indeed include it in our estimate of the 
transverse size of the nucleon in Section 

A comment is in order concerning the overall or- 
der of the pion cloud contribution to the gluon den- 
sity within the \/Nc expansion. If one infers the Nc- 
scaling of g-„NN from the Goldberger-Treiman relation, 
gjrNN = M]\rgA/ Fjr, in which the isovector axial coupling 
of the nucleon, gA, scales as Nc, and the pion decay con- 
stant, Ftt, as Nc^'^ , one is led to g-j^NN ^ Nc^"^ ^3- This 
would mean that the distribution of pions in the nucleon 
at large 5, Ea. l36|) . as a function of y, behaves as 



Nj X function(7Vc?;, b). (44) 



This scaling behavior is analogous to that of the gluon 
distribution in the nucleon in the large-TVc limit . As- 
suming that the gluon distribution in the pion, gjri^z), 
does not scale explicitly with Nc, Ea. H44() would imply 
that the large-iVc scaling of the large-6 contribution to 
the gluon density, as defined by the convolution integral 
(f^ . is of the form 

5(2;, 6)cioud ^ Nc X iunction{NcX, b), (45) 

i.e., the same as for the gluon distribution at average val- 
ues of b. This seems natural, as going to large b should 
not change the TVc-scaling behavior of the gluon distri- 
bution [remember that we consider 6 = 0(1) in l/Nc]. 

Our discussion of pion cloud contributions to the nu- 
cleon's parton distributions in the large-iVc limit so far 
pertained to isoscalar distributions — specifically, the 
gluon distribution. To complete the picture, it is instruc- 
tive to consider within the same approach also isovec- 
tor distributions, which are known to be suppressed rel- 
ative to the isoscalar ones in the 1 /Nc expansion (for the 



polarized distributions, the relative order would be re- 
versed) \XM. The simplest example is actually the flavor 
asymmetry of the sea quark distributions in the proton, 
53(2:) = d{x) — u{x). This distribution has extensively 
been studied within the traditional pion cloud model, in 
which no restriction is placed on impact parameters, but 
the the pion loop integrals are cut off by form factors 
at the ttNN and ttA^A vertices, see Refs.,lX] for a re- 
view. We can, of course, also study it in our approach, 
where pion exchange is considered only at impact param- 
eters b ~ 1/M^. The analogue of Ea. H22|l . describing the 
large-5 asymptotics of the ffavor asymmetry, including 
both N and A contributions, is (see e.g. Ref . ) 



<l3.{x,b) 



<l7r{x/y) 



g /ttJV (?/,&) 



2 f-RA{y,b) 



(46) 



Here qTr{z) = c?7r+(-z) = UTr-{z) is the antiquark distri- 
bution in the charged pions, which may safely be ap- 
proximated by the valence quark distribution, qTr{z) « 
{l/2)v^{z) in the conventions of Ref.|23|. The distri- 
butions of pions in the nucleon, /^at and /^a are the 
isoscalar distributions (sum of the distributions of tt"*" , 7r~ 
and 7r"), as defined in Eas. (|24() and H40|l : the isovector 
nature of q3{x, b) is accounted for by the factors 2/3 and 
— 1/3. Note that here the A contribution enters with op- 
posite sign relative to the N one. In the large-A^c limit 
/^A is exactly twice as large as /^at, see Ea. (|^ . and 
the two contributions in Ea. (|46|l cancel exactly. This is 
as it should be: On general grounds, the isovector un- 
polarized parton distributions are suppressed at large Nc 
compared to the isoscalar ones 14], i.e., they scale as 
Nc X function(iVca::) as opposed to Eg . (|45|l . and if the 
leading pion cloud contribution did not cancel it would, 
following the logic of the preceding paragraph, induce a 
contribution with the same iVc-scaling as in the isoscalar 
case, in violation of the general scaling behavior. Thus, 
we see that, also in the isovector case, the large-6 contri- 
butions due to pion exchange respect the basic 7Vc~scaling 
properties of the parton distributions. 



V. TRANSVERSE SIZE OF THE NUCLEON AT 
SMALL X 

We now proceed to investigate the relevance of the 
asymptotic behavior of the impact parameter-dependent 
distribution at large 6, derived in Sections IIIII and IIVI 
for the x-dependence of the two-gluon form factor of the 
nucleon. For simplicity, we concentrate on the overall 
transverse size of the nucleon, (5^), which determines the 
t-slope of the two-gluon form factor at small t. This will 
reveal an interesting phenomenon — an increase of the 
transverse size of the nucleon at small x due to pion cloud 
contributions. 

The average impact parameter squared is defined as 
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an integral over all values of b, cf. Ea. H13() . When com- 
puting this quantity we have to consider the large-6 con- 
tributions, where the gluon distribution is described by 
Eqs.|(22J) and l|26|) . together with contributions from av- 
erage configurations in the nuclcon, which account for 
the bulk of the gluon distribution. The size of these con- 
figurations is determined by the binding of the valence 
quarks in the nucleon. In a schematic "two-component" 
picture we can write 

J (fb b^ [g{x, 6)buik + 6(6 > 60) g{x, 6)cioud] 



= (62)buik + (6") cloud. (47) 

The 6^-integral in the numerator is computed in two sep- 
arate pieces, while in the denominator we have the total 
gluon distribution (bulk plus cloud) for the given value 
of X. The integral over the pion cloud contribution is 
restricted to values b > bo; the cutoff 60 will be specified 
below. To estimate the "bulk" contribution to (6^) we 
should look to the slope of the elastic "two-gluon form 
factor" of the nucleon, cf. Eqs.0 and p4|l . Unfortu- 
nately, this form factor cannot be measured directly in 
experiment. For a rough estimate we can take the form 
factor of a quark operator which does not receive contri- 
butions from the pion cloud, e.g. the axial form factor of 
the nucleon. This gives an estimate of the "bulk" contri- 
bution to the transverse size of 

(6')bulk ~ ^(r-2)axiab (48) 

where the factor 2/3 arises from converting the "three- 
dimensional" axial charge radius into the "two-dimen- 
sional" average of 6^, cf. Ea. ll4|l and after. With the 
experimental value (r^)axiai = 0.46 fm^ this comes to 

(6')buik ~ 0.3 fm^. (49) 

This simple estimate reproduces well the measured t- 
slope of the J/ il) photoproduction cross section at fixed- 
target energies jia |. 

Consider now the x-dependence of the "cloud" contri- 
butions to (6^), which is computed by integrating the 
asymptotic expression (|22ll over 6. We have seen in 
Section IIIII that the characteristic transverse size of the 
pion cloud contribution to the gluon density is of order 
1/Mjv for X > Mtt/Mn, and becomes of order 1/Af^ for 
X <C Mtt/Mn, see Fig. 13 We thus expect an increase 
in (6^)cioud at small x, which should set in at approxi- 
mately X ~ Aly^/Mjy. Assuming that the bulk contribu- 
tion (6^) bulk does not change much over the region of x 
considered, this would result in an increase in the overall 
transverse size of the nucleon, (6^), see Eq. (|T7j) . 

To investigate this effect numerically, we choose the 
GRV parametrization of the LO gluon distribution in 
the nucleon at a scal e = 4GeV^ [2i|, and the 
parametrization of Ref . |23| for the gluon distribution in 




the pion. As shown in Fig.|Sl the rise of the two distribu- 
tions at small x is similar; the ratio of the pion and nu- 
cleon distributions for x < 0.01 is approximately 2/3, as 
would be expected of gluons generated radiatively from 
valence quarks at a low scale. 

Following the discussion in Section IIVI we include 
in the pion cloud contribution to the transverse size 
also contributions from A intermediate states, which, al- 
though theoretically subleading at large 6, are expected 
to be important quantitatively. We define the total dis- 
tribution of pions in the nucleon as in Eq. H39() , and eval- 
uate /ttat and /^a using the asymptotic expressions at 
large 6, Eas. l(^ and (^UJ). We use the physical values 
for the N and A masses; for the coupling constants we 
take gTrNN — 13.5 and gTiNA — (3/2) gTrjvw, which is 
close to the phenomenological value 30]. We then com- 
pute the 6^-weighted integral of the distributions with a 
lower cutoff bo, cf. Ea. (|47|l . whose value we take to be the 
size of the bulk of the gluon distribution in the nucleon, 
60 = (^^)buik ~ 0.55 fm. The integrated distributions 

J dHb" f^N^.A{y,b), (50) 

b>bo 

are shown in Fig. (in units of 69), as functions of y. 
One sees that the maximum of the N contribution is 
approximately at y ~ Mt^/M^. The A contribution sets 
in at smaller values of y compared to — a consequence 
of the different y-dependence of the exponential decay 
constants, Eas. H26() and H41|) — but eventually becomes 
larger than the A^ contribution at smaller y, in qualitative 
agreement with the large- A^c relation, Ea. (|43|l . 

The results for (6^)cioud are presented in Fig. [T] The 
dashed line shows the result for (6^)cioud obtained from A^ 
intermediate states only, cf. Ea. (|24|l . The solid line shows 
the sum of the contributions from N and A intermediate 
states (with physical A^ and A masses), cf. Ea. (|40|l . One 
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FIG. 6: The 6^-weighted integral of the distributions of pions 
in the nucleon, Eg. 15011 . with a lower cutoff 6o = (''^)buik ~ 
0.55 fm, as a function of the pion momentum fraction, y. The 
values of the integrals are given in units of &o- Shown are the 
contributions from (solid line) and A intermediate states. 
These integrals enter in the calculation of the pion cloud con- 
tribution to the transverse size of the nucleon, (&^)cioud, c/. 
Eg.lHTl. 



FIG. 7: The pion cloud contribution to the average transverse 
size squared of the gluon distribution in the nucleon, (6^)cioud, 
Eg. (1471 . as a function of x. The values shown here were ob- 
tained by integrating over the large-6 tail of the gluon distri- 
bution, Eg. 1221 . with a lower cutoff bo = {b'^)\^-^n^ ~ 0.55 fm. 
Dashed line: Contribution from A'' intermediates states, c/. 
Eg. 1241 . Solid line: Sum of A*' and A contributions, c/. 
Eg. 1401 . These results should be compared with the "bulk" 
contribution to the nucleon size, {6^)buik ~ 0.3 fm^. 



sees that the onset of the A contribution is "postponed" 
to smaller values of x compared to the N one, and that at 
X ~ 10~^ the A contribution makes about half of the to- 
tal result These results for (6^)cioud shown in Fig.[7| 
should be compared with the bulk contribution to the 
nucleon size, (6^)buik, which we estimated at ~ 0.3 fm^, 
c/. Eas. H48|l and (|49|l . We see that, with our estimates 
of the parameters, the increase in the total transverse 
size of the nucleon due to pion cloud contributions for 
X <C Mt^/Mm amounts to ~ 20% of the total transverse 
size of the nucleon. This is in qualitative agreement with 
the analysis of Ref. 0|, which concluded that the nu- 
cleon two-gluon form factor is similar to the nucleon axial 
form factor for 0.3 > x> 0.05, and which suggested that 
a large part of the observed increase of the t-slope of the 
J/^ photoproduction cross section between fixed-target 
energies [E^^ < lOOGeV, a: > 0.05) and HERA energies 
{x < 10~^) is due to the contribution of the pion cloud. 

In the results shown in Figs. |H| and [7| the distributions 
of pions in the nucleon f^^N and /^a were evaluated using 
the leading asymptotic expressions at large b, Eas. (|24|) 
and (|40() . which contain only the leading power of 1/(k&) 
in the pre-exponential factor. In principle, since we are 
integrating the distributions over b from a lower cutoff 

^0 — (^^)buik ~ 0.55 fm, which is numerically compara- 
ble to (in fact, even smaller than) 1/(2M^) — 0.73 fm, 
subleading terms of order 1/(k6) in the pre-exponential 
factors, which become ~ 1/(2M^6) at small y, could be- 
come important. To estimate the influence of these terms 
we have evaluated the integrals (|50|l also with the full 
expression for f^^N derived from the pion exchange dia- 
grams, Ea. ljA19|) . and the corresponding result for /tta 



obtained by numerical evaluation of Ea. l|A29() . Wc find 
the 6^-weighted integral (|Sn|l in this case to be ~ 30% 
larger for the N, and ^ 10% larger for the A contribu- 
tion, resulting in a 20% larger total value of (&^)cioud 
at X <C Mt^/Mm than the one shown in Fig. [7| (solid 
line). The small difference compared to the results ob- 
tained with the leading asymptotic expressions for the 
distributions at large b clearly demonstrates that the ef- 
fect we are considering is indeed due to large transverse 
distances, where our approximations are justified. Note 
that this is because we are considering integrals over the 
distribution of pions weighted with b^, which strongly en- 
hances contributions from large distances; it would not 
be the case for the simple integral of the distribution of 
pions over b (the number of pions with b > bo) which 
would enter in the pion cloud contribution to the parton 
density, such as the flavor asymmetry of the sea quark 
distributions, cf. the discussion in Section Hvl To con- 
clude, we emphasize that the uncertainty discussed in 
this paragraph concerns only the height of the "jump" 
of (6^)cioud between x > Mj,/Mn and x < M„/Mn; 
the basic mechanism of suppression of the pion cloud 
at a; > Mt^/M^ is due to the exponential factor in the 
asymptotic expansion at large 6, and thus more robust. 

In the large-iVc limit, the sum of the N and A con- 
tributions to (&^)cioud would be exactly 3 times larger 
than the contribution alone. As one sees from Fig. [7| 
our result for the sum of N and A contributions with 
physical N and A masses (solid line) lies below 3 times 
N only (dashed line) by a factor of ~ 2/3. It thus ap- 
pears as if the large-A^c limit overestimated the "exact" 
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result. However, there is the question which value one 
should take for the common N and A mass in the large- 
Nc limit. Multiplying the N contribution in Fig. [3 by 
3 amounts to taking Ma — Mm — 940 Mev. This is an 
extreme choice, and the resulting (6^)cioud should clearly 
be seen as an upper limit. Choosing larger values for the 
the common mass in the large- A^c limit, as are suggested 
by most chiral soliton models, would lower the large- A'c 
estimate considerably, bringing it more in line with our 
"exact" result. 

In the derivation of the large-6 asymptotics of the 
gluon density in the nucleon, Ea. H22() . we neglected the 
"intrinsic" transverse size of the gluon distribution in the 
pion. Of course, the transverse size of the distribution in 
the pion is expected to grow at small z = x/y (momen- 
tum fraction of the gluon relative to the pion). This does 
not play any role in the region of the onset of the growth 
of the transverse size of the nucleon, x ^ M^^/Mn^ since 
the configurations responsible for the growth of the nu- 
cleon size are those at the lower limit of the y-integral, 
y ^ X for which z ~ 1. However, it becomes an issue 
for values of x considerably smaller than M-^/Mpf. In 
other words, the finite size of the pion does not change 
the onset of the growth of the nucleon size in x, but may 
increase the height of the jump. 

One could think of applying the reasoning of Section lHll 
to the pion itself and consider the growth of the trans- 
verse size of the pion due to soft-pion exchange in the 
i-channel. In the case of the pion the growth of the 
transverse size should start already for values z ^ 1, 
contrary to the nucleon, where it is "postponed" until 
X ^ Mj^/Mn. However, since the pion-pion coupling is 
weaker than the pion-nucleon coupling the relevance of 
the "pion cloud of the pion" for the transverse size of its 
gluon distribution is questionable. In particular, in the 
large-A'c limit the pion-pion interaction is suppressed 
relative to the pion-nucleon interaction. The transverse 
size of the pion is thus likely to be determined by its 
quark core and can be studied only in models. A rough 
estimate of the transverse size of the pion is provided by 
its electromagnetic charge radius [cf. Eq.l^SJ] 

- ^{r^). = 0.3 fm^. (51) 

This suggests that the size of the gluon distribution in 
the pion is comparable to that of the bulk of the gluon 
distribution in the nucleon. The finite size of the pion 
gives a positive contribution to (6^)cioud7 which in prin- 
ciple should be added. However, this would exceed the 
accuracy of our simple two-scale picture in which the core 
radii (either nucleon or pion) are parametrically smaller 
than 1/M^. A quantitative estimate of the effects of 
the finite transverse size of the pion can be made within 
the chiral quark-soliton model of the nucleon, which de- 
scribes the structure of the pion and the nucleon in a 
unified way [T3 . [T5l |. This approach would also incor- 
porate the effects of the pion-nucleon form factors in a 
consistent way. 



VI. PION KNOCKOUT IN HARD EXCLUSIVE 
PROCESSES 

We now turn to the question whether the pion cloud 
contributions to the gluon density in the nucleon at large 
impact parameters could directly be observed in exclu- 
sive photo/electroproduction experiments, and whether 
in this way one could extract information about the gluon 
distribution in the pion. 

One possibility would be to measure the i-dependence 
of the cross section for p electroproduction off the nu- 
cleon, 7*-|-A^ — > p+N, or for J/tp~ or T photoproduction, 
over a sufhciently wide range such that one can restore 
the 6-dependence for b > 1/(2M^) » 0.7 fm. (A simi- 
lar program was carried out in Ref . |3ll | in the framework 
of the dipole picture of high-energy scattering.) In this 
way one could extract the gluon density in the nucleon at 
large 6, which could be compared with the asymptotic ex- 
pression Eas. H22|) . H24|l and H2t)|) . A problem with this ap- 
proach is that the theoretical prediction Eq. (|22|l has the 
form of a convolution in the momentum fraction of the 
pion, y, which would render the extraction of the gluon 
distribution in the pion difficult. Another problem is the 
unknown size of sub-asymptotic contributions to the nu- 
cleon's gluon distribution at large b, which have a faster 
exponential fall-off than the N contribution but may be 
enhanced by numerical prefactors. On the practical end, 
the resolution in t at small t, and the t-range covered 
by the present HERA experiments, are definitely insuffi- 
cient for making the transformation to 5-space with the 
required accuracy. Hopefully, the necessary resolution 
could be reached at the planned Electron-Ion Collider 
(EIC) 32]. 

A more direct measurement of the glue in the pion 
cloud is possible in exclusive reactions in which the 
"struck" pion is observed in the final state. Here we wish 
to consider processes of the type 

J* + N a + TT + N' (or A), (52) 

in the kinematical region where the subprocess 7* + tt — s- 
a + TT satisfies factorization theorems for exclusive pro- 
cesses, see Fig. |SJi. This could be either a large-Q^ pro- 
cess, with the produced particle a a real photon or a 
vector meson (in the latter case only the amplitude for 
longitudinal polarization of the virtual photon should be 
considered), or a process in which a heavy quarkonium 
is produced (a = J/'0, '^)- The process (I52II is an exam- 
ple of a general class of reactions in which only a "white 
cluster" of configurations in the target nucleon actually 
participate in the hard scattering. It was pointed out 
in Ref. [s^ that such processes can be used to obtain 
qualitatively new information about baryon structure. 

To describe the process (I52f) we choose p and q collinear 
and along the 3-axis. The outgoing particles are then 
characterized by their longitudinal momentum fractions 
and transverse momenta. The momentum of the ex- 
changed pion is completely fixed by the external mo- 
menta and thus measurable; in particular, fc^ = (p'— p)^. 
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FIG. 8: (a) The mechanism for hard exclusive photo- or elec- 
troproduction of a system a (vector meson, heavy quarko- 
nium), accompanied by wide-angle pion production, — + 
a+n+N' (or A). The blob denotes the amplitude for the sub- 
process 7 + TT ^ a + TT. (b) Competing reaction mechanisms 
in which the whole nucleon participates in the hard process, 
and the pion is emitted from the external baryon lines. These 
contributions can be suppressed if the pion is emitted under 
sufficiently large transverse momentum. 



We characterize the momentum of the exchanged pion by 
its longitudinal momentum fraction y and transverse mo- 
mentum kx • The ampUtude for the process l|52|) is given 
by the product of the the amplitude to find the pion in 
the nucleon with momentum fraction y and transverse 
momentum = — p^, and the amplitude of the pro- 
cess ^ + TT ^ a + TT. At sufficiently small Bjorken x the 
latter process is dominated by the gluon distribution in 
the pion. The ^-dependence of the latter is given by the 
two-gluon form factor of the pion, TT^(t). As we argued 
already in Section HTTI the characteristic scale of this form 
factor, ni2g, is comparable to that of the pion electromag- 
netic form factor, and the analogy with the nucleon 
suggests that it may be as large as m^^ « 1 GeV^ , see 
Ea. (|17|l and after. In any event, the characteristic scale 
of this form factor is much harder than the transverse 
momenta which could be generated in the pion-baryon 
cluster system. This is crucial for our discussion, since 
it guarantees that the reaction mechanism of interest, 
Fig. can be separated from the competing mecha- 
nism in which the whole nucleon is involved in the hard 
process, and the extra pion is emitted by the incoming 



or outgoing nucleon, see Fig.|S|3. Keeping {p — p'Y small, 
of order Af^, and choosing t — {V — g)^ sufficiently large 
(|1±| > 0.3 . . . 0.4 GeV), one effectively suppresses contri- 
butions of the type of Fig. ^p. 

We can write the the modulus squared of the amplitude 
of the process (|52f) in the form 



|M 



■y+N^a+TT+N' I 



X 



47r 
V 

\M. 12 



7-|-7r^a+7r | 



(53) 



Here wn^n' is an isospin factor; its value is 1 for p ^ p 
or n — > n, and 2 for p n or n p transitions. Fur- 



thermore, Ai^+TT^a+TT denotes the invariant amplitude of 
the subprocess j + tt a + it. The function d-^NiUi k_L) 
is defined as, up to a factor, the modulus squared of the 
amplitude for the nucleon to emit a pion, including the 
pion propagator, averaged (summed) over the incoming 
(outgoing) nucleon helicities: 



47r 



rf7r7v(2/, kj 



1 91nn \u'i-l5T°-u\'^ 

2 ^ (fc2 - M2)2 ■ 



Explicit calculation gives 

d^Ar(j/,k^) = 



2 

SttNN 



(54) 



(55) 



where s = {y^M^ + k^)/y = — fc^ ig the spacelike 
pion virtuality. This function can be interpreted as the 
transverse momentum-dependent distribution of pions in 
the nucleon. Formally, the integral of Ea. H55|) over all 
k^, multiplied by 3 to account for the isospin degener- 
acy, gives the total distribution of pions in the nucleon, 
H^iy,-A'i ^ 0), cf. Eq.JSnil for A_l = 0. Again, it is 
physically more sensible to consider the distribution of 
pions H55|l in coordinate space. We define the Fourier 
transform of Eq. H55(l as 



d^Ar(?/,a_L) = 



e'(''^"^) d.iv(y,k^), (56) 



where is a transverse coordinate variable conjugate to 
k^. At large |a^| it behaves as 



-KN\a±\ 



where 



^ym% + yMl 



(57) 



(58) 



For y ~ Mt^ /Mn the width of this distribution in trans- 
verse space is of the order 1/Af^, i.e., much larger than 
the size of typical configurations in the nucleon, 1/Mjv. 
In this region the exponential dependence could in prin- 
ciple be observed experimentally in the Fourier trans- 
form of the cross section of the reaction H52|) with re- 
spect to k^ = — p'j^, i.e., the transverse momentum of 
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the outgoing nucleon. Due to the broad distribution of 
the cross section in t it should be easy to make a cut on 
—imin which would allow to measure dj^Niv, k±) in a suffi- 
ciently wide range of p' j_ such that one could reconstruct 
d^Ar(j/, a_L). Alternatively, one could model d^Ar(y,aj^) 
in the whole range of |aj_|, calculate the Fourier trans- 
form, and check how small needs to be in order for 
the region of |a^| which we consider as safe to give 
the dominant contribution. We assume here that for 
|k^| ^ jk'^l the dependence of the amplitude of the sub- 
process 7-|-7r— > a-|-7ron the transverse momentum of 
the "initial" pion, k^, can be neglected. 

At sufficiently small Bjorken x the hard subprocess 
7 + 7r^a + 7ris dominated by the gluon distribution 
in the pion. Measuring the cross section of the process 
7-I-7V— > a-|-7r-|-iV', and isolating the contribution from 
pion exchange as described above, one could thus mea- 
sure the gluon distribution in the pion, including its t- 
dependence (two-gluon form factor). The complete dis- 
cussion of the cross section of this process, appropriate 
kinematical cuts, etc. is beyond the scope of the present 
investigation and will given elsewhere. Here we note only 
that at sufficiently small x a simple expression can be 
written for the ratio of the squared amplitudes of the 
process and the corresponding production process 
without pion emission: 



\M. 



47r 

y 



y+N^a+-iT+N'\ J i' 1 N 

— = WN^N' d^N[y,iiA.) 



xg{x, Q^) 



n 2 



(59) 



Here denotes the scale appropriate for the hard pro- 
cess, see Refs.0,13 for a discussion, and TTj{t) and T{t) 
are the two-gluon form factors of the pion and nucleon, 
respectively (in both processes t denotes the invariant 
momentum transfer between the incoming photon and 
the outgoing vector meson/photon, see Fig. |Sli). Since 
the distribution of pions is not too small for small y, and 
the gluon density in the pion strongly increases with de- 
creasing x/y, we expect the cross section for the process 
(|52|l to be non-negligible for y ^ 0.05 ... 0.1. In addition, 
if indeed the two gluon form factor of the pion decreased 
slowly with t as in Eq. (|17|l . the production rate would 
decrease with jk'^j only as cx l/jk'^p. One may thus ex- 
pect noticeable rates for pions produced with transverse 
momenta |k^| > 2 GeV, which is within the acceptance 
of the current HERA detectors. 



VII. CONCLUSIONS AND OUTLOOK 

In this paper we have studied the role of the pion 
cloud of the nucleon for the gluon density at large trans- 
verse distances. We have identified a specific contri- 



bution to the gluon density of transverse size 1/M^, 
due to soft-pion exchange, which becomes "visible" for 
X ^ Mtt/Mn. This results in a finite increase of the aver- 
age transverse size of the nucleon as x drops significantly 
below M-^/Mm, in agreement with the observed change 
of the i-dependence of the J/i/; photoproduction cross 
section with incident energy. While this phenomenon 
has been described previously at the qualitative level, the 
present investigation shows that with the phenomenolog- 
ical parameters for soft-pion exchange one indeed obtains 
a consistent quantitative picture. 

Let us summarize again in which range of x the soft- 
pion exchange contribution to the gluon density is mean- 
ingful. We have seen that this mechanism sets in for 
X < M^/Mjv ~ 10~^, when the relativistic kinematics 
allows the pions to propagate distances of order l/M^ 
in the transverse plane. It grows to its full strength for 
X ^ 10~^, when the momentum fraction of the gluons rel- 
ative to the pions becomes small (~ 10~^). In principle 
this contribution to the gluon density survives down to 
considerably smaller values of x. However, for such val- 
ues of X other effects such as multi-step Gribov diffusion 
become dominant. 

The gluon density at large transverse distances is an 
important element of the structure of the nucleon which 
can be calculated theoretically in a model-independent 
fashion. It can be probed in exclusive processes such as 
J/ip photoproduction, or hard electroproduction of real 
photons or vector mesons at small x, by isolating con- 
tributions from large impact parameters through mea- 
surement of the ^-dependence. The QCD factorization 
theorem guarantees that this feature of the nucleon is 
process-independent as long as the conditions for factor- 
ization arc met. A particularly promising way of access- 
ing the glue at large transverse distances appears to be 
hard exclusive production off the nucleon with associated 
"pion knockout" , which should be possible to observe at 
HERA. 

Our results for the asymptotic behavior of the gluon 
density at large transverse distances provide model- 
independent constraints for parametrizations or dynami- 
cal models of the impact parameter dependence of the nu- 
cleon's gluon density. Such parametrizations are needed 
to describe e.g. diffractive Higgs production in proton- 
proton collisions at LHC energies [23| ■ 

The techniques developed in the present paper can 
readily be extended to quark distributions, including po- 
larized distributions. The large-6 limit of these distri- 
butions is in principle also due to soft-pion exchange, 
and can be described in a model-independent way. How- 
ever, the minimum number of exchanged pions in the 
i-channel depends on the G-parity of the operator mea- 
suring the parton distribution, and the relevance of the 
asymptotic region for eventual observable physical phe- 
nomena needs to be studied case-by-case. The impact 
parameter representation of parton distributions is also 
useful in connection with specific dynamical models of 
the nucleon, as it leads to new insights into the dynam- 
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ical origin of these quantities. For example, with regard 
to the flavor asymmetry in the nucleon's sea quark dis- 
tributions, u{x) — d{x) and Au{x) — Ad{x), it may help 
to resolve the long-standing issue to which extent these 
asymmetries are due to the pion cloud of the nucleon, or 
to Pauli blocking by the valence quarks (36l |. 

In the present investigation we have limited ourselves 
to the generalized parton distributions at zero longitu- 
dinal momentum transfer (zero "skewedness" ) , which is 
sufficient for studying the impact parameter dependence 
of the usual parton densities in the nucleon. The inclu- 
sion of skewedness effects is an important problem which 
we leave for future treatment. In particular, this would 
allow to include into the considerations a much wider 
range of experiments, such as deeply- virtual Compton 
scattering or meson production (pseudoscalar and vec- 



tor) in the valence region. First steps have been taken 
to extend the impact parameter representation to gener- 
alized parton distributions at non-zero skewedness [35l |: 
however, one is still lacking a simple interpretation as is 
available in the diagonal case. 
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APPENDIX A: DISTRIBUTION OF PIONS IN 
THE NUCLEON AT LARGE IMPACT 
PARAMETERS 

1. Evaluation of the pion exchange graphs 

In this appendix we outline the derivation of the 
asymptotic behavior of the "parton distribution" of pi- 
ous in the uucleon at large impact parameters. We first 
consider the distributions obtained with N intermediate 
states only; the A contributions will be included subse- 
quently. 

In the limit of large transverse distances it is justified 
to treat the nucleon and the pion as pointlike elemen- 
tary particles. We assume a pion-nucleon coupling of 
pseudoscalar form 



NN 



(Al) 



where r" are the isospin Pauli matrices. The result for 
the large-6 behavior is in fact independent of the type of 
coupling; i.e. it is the same for all couplings which are 
equivalent in the soft-pion limit. 

We first evaluate the ^-dependent distributions of pi- 
ons in the nucleon, as defined by Ea. (|21|l . The nucleon 
matrix element of the pion light-ray operator is com- 
puted using standard Feynman rules, 

{N{P + A/2)| 7r'^(-An/2)7r'^(An/2) \N{P - A/2)) 



o ■ 2. 



NN 



(2^ 



^—iX{kn) I iX(kn) 



[{k + A/2)2 - M2 + lO] [{k - A/2)2 -M^ + iO] 
u' [P — k + Mpf) u 



[{P - fc)2 - Af 2, + lO] 



(A2) 



The factor 3 arises from the summation over the pion 
isospin projection. In order to read off the expressions for 
the distribution functions we have to convert the Dirac 
structures in the numerator to the form of the L.H.S. of 
Ea. (|21() . The 75 matrices can be eliminated using the 
anticommutation relations. The integral over the k term 
can be projected on the structures P, A and n by appro- 
priate projections of the four-vector k. Finally, we use 
the identities u'Pu = Mmu'u and u' Ku = 0, as well as 



3l 

M. 



1 



N 



U 'Ji.U - 



2Mn 



(A3) 



which follow from the Dirac equation for the external 
nucleon spinors. The net result is that wc can replace in 
Eq.lEU 



u i"f5 (P 
1 

(P^ 
2(Pn)2 



(fcP) 



Mn ) hb u 

+ (p2 



(kn) 
(Pn) 



u nu 



(A4) 
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To get the distribution of pions as a function of the mo- 
mentum fraction y we still have to compute the Fourier 
transform of Eq. l)A2p with respect to the longitudinal dis- 
tance variable A, see Ea. (|21|l . The integral of the com- 
bined exponential factors results in delta functions, 



2tt 



S[y{Pn)T(kn)]. (A5) 



Combining everything, we see that the distribution of 
pions resulting from TV intermediate states is given in 
covariant form by 



H^Niy,t) = 6iglf^^ {Pn)y 



_cffc 



j6[y{Pn) - (kn)] 





A/2)2 - 


M2 H 






1 






[{k- 


A/2)2 - 








1 







[(P - fc)2 -M^ + lO] 
(kn) A2 



ikP) - 

{y -y), 



(Pn) 4 



(A6) 



where t — A^. For the helicity-flip distribution E^rNiy, t) 
the expression in brackets in the next-to-last last line 
should be replaced by M"^ (kn) / (Pn) . In the following we 
quote only the results for the distribution i/^Ar used in 
the present investigation; the corresponding expressions 
for E-^N can be derived analogously. 

To evaluate the momentum integral we choose the 
frame described in Section^ where — (1,0,0,-1), 
nP = P+, and = 0, sec Eqs.Q, which implies 
P- = P^/P+. Furthermore, A+ = A" = 0, and thus 
A'^ = —t. Corresponding light-like components are in- 
troduced for the loop momentum, k. The integral over fc"*" 
is taken trivially; the delta function restricts fc+ to yP~^. 
The integral over k~ is then performed using Cauchy's 
theorem. The poles of the three propagators are located 
at 



(k^±Ax/2)2 + A42~zO 
yP+ 

-kj - yAf^ + yAl/A + ^0 
yP+ 



(A7) 
(AS) 



The integral is non-zero only if the poles lie on different 
sides of the real axis, which implies y?; > 0, or < y < 1. 
Closing the contour in the upper half-plane, encircling 
the pole ljA8|l . we obtain for < y < 1 



H^N{y,-Al) = 



^91nn y 



»7r 



y J (2^)2 

s_ - yA\ 



(A9) 



where 



s± 



[ik^±yA^/2f + y^M^^] 



-{k ± A/2f at the pole (AlO) 



are the pion virtualities. The minimum value of the pion 
virtuality in the loop is 



y'^M] 



N 



y 



(All) 



For A^ = our result (jA9|l reduces to the well-known 
expression for the isoscalar distribution of pions in the 
meson cloud model, f-T^Niy) of Ref.fl3|, up to form factors 
at the pion-nucleon vertices, which are not needed in our 
approach, see below. Note that in Ref.^3 the isospin 
factor 3 is included in the definition of Qt^nn- 

In Ea. ljA9|) it is assumed that < y < 1; the 
distribution for — 1 < y < follows trivially from 
H-^wi-y, -^D = -H^N{y, -A5_). In the following, for 
simplicity, we shall always assume that < y < 1, and 
that the distribution at negative y is restored in this way. 

The transverse momentum integral in Ea. ljAQI) for the 
distribution HT^{y,—A\) is logarithmically divergent. 
However, this divergence is irrelevant for the large-6 be- 
havior of the corresponding impact parameter-dependent 
distribution. It could be removed by subtracting the 
function at — A^ = i = 0, 



H^N{y,-Al_) 



H^N{y,-Al)-H^N{y,Q). (A12) 



In impact parameter space this would amount to sub- 
tracting a delta function at b = 0, which is irrelevant at 
large h. With the saddle point method employed below, 
we do not need to perform this subtraction explicitly. 

For deriving the large-6 behavior of the correspond- 
ing impact parameter-dependent distributions it is con- 
venient to introduce the integral representation 



1 



da e 



-otA 



(A>0) 



(A13) 



for the kj^ -dependent denominators in Ea. l|A9|l . The in- 
tegral over kj_ then reduces to a Gaussian integral, which 
can be computed. We obtain 



i?^Ar(y,-Ai) 
_ 



167r 
X exp 
1 



da I df3 

"'0 



a + P 



-(a + ^)(Smin,W + M2) 



^min,A'' 



a/3 



a + P 



a + P 

-A 2 

y^i 



{a + PY 



,(A14) 



where Smin.A' is defined in Ea. (|Alip . We can now pass 
to the impact parameter-dependent distribution of pi- 
ons in the nucleon, c/. Eq. H23|) . The Fourier transform of 
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Ea. (jA14|l with respect to Aj^ is again a Gaussian inte- 
gral. We get 



UN{y,h) 

__ "^drrNN y_ 

Q47J-3 y 



da 



dp 



X exp -{a + (3){s 

+ Ml) 

If 1 62\ 

ap \ ap Ay J 



aP Ay 



(A15) 



In this representation the large-6 asymptotics can be 
studied using the saddle point method. The stationary 
point of the exponent is at 



= P - ^[H^ 



..N 



Mi 



-1/2 



(A16) 



and one immediately sees that with exponential accuracy 
the large-& asymptotics is given by 



-KNb 



with 



KN 



Ml 



V 



1/2 



(A17) 



(A18) 



One can easily determine the pre-exponential factor by 
computing the Gaussian integral over deviations from the 
stationary point; the result is given in Ea. H24|l . Alterna- 
tively, one can express the integral Eg. ljAlSjl in terms of 
modified Bessel functions, 



f7vNiy,b) 



^91nn y 



|smiiKiV [KQ{KNh/2)Y 



IGtt'^ y 

+ (s,ni„,jv + Af2) [ifi(Kw6/2)]'} , (A19) 



and recover the asymptotic behavior Ea. (|24|l by substi- 
tuting the asymptotic expressions for the Bessel functions 
at large arguments. In particular, in this way one can get 
also the subleading corrections to the pre-exponcntial fac- 
tor involving higher powers of l/{K]\[b), which arc more 
difficult to compute using the saddle point method. 



2. Large— & behavior from singularities 

The large-& behavior of the distribution of pions in the 
nucleon can also be derived in an alternative way, which 
directly exhibits the connection of the large-6 asymp- 
totics with the singularities of the pion and nucleon prop- 
agators in the Feynman integral (jA6(l . Performing in 
Ea. (|23|l the integral over the angle between the trans- 
verse vectors and b we obtain (here = |A^|) 



UNiy,b) = 



rfAi 
~27 



Mb A J 



i?,Ar(y,-Ai), 
(A20) 



where Jq denotes the Bessel function. Substituting 
the latter by its asymptotic form for large arguments, 
MbA±) - (27r6Ai)-i/2 [e*(6A^-^/4) ^g-»(bA^-^/4)]^ 

and rotating the integration contour in the complex A^ 
plane to run along the positive (negative) imaginary axis 
in the integral with the first (second) exponential factor, 
setting = ±iS + with S > real, we get 



f7vNiy,b) 



1 



dsVs, 



-bS 



threshold 

X -lmHT,N{y,S^ + iO). 

TT 



(A21) 



The large-& behavior of the impact parameter-dependent 
distribution is governed by the threshold behavior of the 
imaginary part of the ^-dependent distribution of pions, 
HTrN{y,t), which appears for some nonzero positive t. 
The imaginary part can directly be computed from the 
Feynman integral ljA6|) using the Cutkosky rules, replac- 
ing the propagators by appropriate delta functions. One 
finds that the i-channel cut starts at d = kn, and that 
the leading threshold behavior of the imaginary part for 
(5 ^ ktv + is 



- lv[\H^N{y,5'^ 

IT 



iO) 



-ig o - (2Sinin,Ar + M^j 

Idtt^ y 



_K%{S - kn) 



1/2 



(A22) 



Substituting this in Eq. (|A21|I and computing the integral 
over S, one obtains the asymptotic behavior quoted in 



3. The A contribution 

The calculation of the contribution from A intermedi- 
ate states proceeds along the same lines as that of the 
N contribution; we only outline the essential steps. We 
assume a ttA^A interaction Lagrangian of the form 



Ltt- 



NA — 



V2M 



N 



3^9^^° A+ 



p a^7r+ A° 
V3 ^ 



(A23) 



where p is the proton field (we skip the terms for the 
neutron), and A^ the spin-3/2 Rarita-Schwinger field 
describing the A resonance [not to be confused with the 
four-vector of the momentum transfer in the matrix el- 
ement (|A2ll ]. Our definition of the coupling constant 
g-rrNA corrcsponds to that of R.ef. . and differs from 
the one used by Koepf et al. [3, (sf,r]VA)Koepf ot ah = 
\/2 (.g:,TAfA)this work- The relative coefficients of the terms 
in Ea. (|A23|l are dictated by isospin invariance. One sees 
that the total distribution of pions in the proton (i.e., the 
sum of distributions of 7r"'',7r~ and tt'^) resulting from A 
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intermediate states is equal to the tt A"*"*" contribution 
multiplied by (1 + 2/3 + 1/3) = 2. 

It the calculation of the Feynman graphs for the A con- 
tribution, we use the propagator of the spin-3/2 Rarita- 
Schwinger field in the form (I is the momentum) 



I + Ma 



Ml 



1 



3Ma 



(A24) 



After projecting on the Dirac structures appearing in the 
L.H.S. of Ea. (|21() . in analogy to the N case, we get [c/. 
Eq.dXSJ] 



9^ [Pn)V 

iM%Ml ^ 

1 



(2^ 



5[v{Pn) - {kn)] 





A/2)2 - 
1 








A/2)2 - 
1 


M2H 


hiO] 



[{P - kf - Ml + iO] 
X [{P ~ k, kf - (fcA)V4 + Ml{~e + AV4)] 
[kn) A2 



-{kP) - 



(Pn) 



Mn{Mn + Ma) 



(y -y)- 



(A25) 



For the helicity-flip distribution -BttaIj/iO the expres- 
sion in the bracket in the next~to~last line should be 
replaced by [Mfjikn)/ (Pn) - MNiMN + Ma)] ■ Intro- 
ducing light-like coordinates as above, and taking the 
integral over k~ using Cauchy's theorem, we finally get, 
for < y < 1 [cf. Eq.J^] 



if^A(y,-Ai) 



o2 



d-k 



ASnM^MlyJ (2^)2 (s+ + M2)(s_ + M^) 
[2Mi(,5+ + ,s_- Ai) 



+ (s+-Afi + A4)(5_ 



Ml 



- yAl] , 



X [2(AfA + AfAr)^ + S. 

where now the pion virtualities are given by 



s± = i(ki±yAx/2)2 

y 



(A26) 



(A27) 



with 



^min, A 



y{Ml - Af^)] 



(A28) 



For zero momentum transfer, Aj^ = 0, Ea. (jA26p reduces 
to the A contribution to the isoscalar distribution of 
pions in the meson cloud model quoted in Ref . ^1_8] [cf. 
the comment on the different conventions for ^tttva after 
Ea. (|A23|l ]. up to form factors at the ttA^A vertices, which 
are not needed in our approach. Again, we are interested 
in Eq. ljA26(l only as a means to derive the large-6 asymp- 
totics of the corresponding impact parameter-dependent 
distribution. This is done in complete analogy to the N 
case, using the integral representation (|A13|I to convert 
Ea. (|A26|l to a Gaussian integral. For reference, we quote 
the result for the 6-dependent distribution, correspond- 
ing to Ea. (|A15|l : 



UAiy,b) 



y-irNA 



iM^m^Ml y 7„ io '^^ 



-{a + P){Sn,in.A+M^)- 



X exp 



2, a + P b"^ 



{X + Al){X + Al) 



a(3 



where 



X 



XA+A_ 



9 



(A29) 



Al 



Aal3y 

,A + {Ma ± Afjv)' 



'-'mill 



^ (A/a ± vMn)^ 

y 



(A30) 



(A31) 



The ellipsis in Ea. ljA29|) denotes terms which produce 
only subleading corrections in 1/6 to the pre-exponential 
factor in the large-6 asymptotics, which we have not writ- 
ten out for brevity. The large-6 behavior can again be 
studied using the saddle point method; the exponential 
decay constant is given by Ea. ljA18|) with Smin.N replaced 
by Sniin.A- The result for the leading large-6 behavior 
is quoted in Ea. H40|l . Note also that the representa- 
tion l|A29l) can be used for direct numerical evaluation 
oi UA{y,b). 



